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Abstract 

We suggest predictive scenarios for neutrino masses which provide a common origin for 
CP violation in early universe cosmology and in neutrino oscillations. Our setup is the seesaw 
mechanism in the context of MSSM with two quasi-degenerate right-handed neutrinos, with 
baryon asymmetry generated via resonant leptogenesis. Three different models are found 
with specific textures in the Yukawa coupling matrices, each with a single phase which controls 
leptogenesis and neutrino CP violation. One model leads to normal hierarchy of light neutrino 
masses and the prediction tan ^13 ~ sin 612 resulting in a value of the reactor mixing angle 

^13 very close to the current experimental lower limit. The other two models predict inverted 
hierarchical neutrino mass spectrum with the sum rules sin^ 612 — ^ — tan 623 sin ^13 cos 6 and 
sin^ 612 — ^ + cot ^23 sin ^13 cos 6 respectively. We obtain a lower bound for the phase \6\ in the 
normal hierarchical model, and a narrow range for \5\ for the inverted hierarchical model from 
cosmology. In our scenario, the mass-splitting between the quasi-degenerate right-handed 
neutrinos arise via renormalization group flow, which provides a lower limit on the MSSM 
parameter tan/3 > 12. The right-handed neutrino masses can be as low as TeV, which would 
avoid the gravitino problem generic to supersymmetric models. 
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1 Introduction 



While the standard model (SM) of strong and electroweak interactions has been extremely successful 
in confronting experimental data, it leaves several questions unanswered. On the observational 
side, it does not provide a viable dark matter candidate, nor a dynamical mechanism to explain 
the observed baryon excess in the universe. Furthermore, the model needs to be extended, albeit 
in a minor way, to accommodate small neutrino masses as needed for atmospheric [1] and solar 
neutrino oscillation data [2] . On the theoretical side, the model suffers from the quadratic divergence 
problem, which destabilizes the Higgs boson mass. 

An elegant synthesis of these issues is provided by low energy supersymmetry (SUSY) and the 
seesaw mechanism [3]. Low energy SUSY can cure the quadratic divergence problem for the Higgs 
boson mass. In its simplest form it also provides a natural dark matter candidate, the lightest SUSY 
particle (LSP). The seesaw mechanism assumes the existence of right-handed neutrinos (RHN) 
N which facilitates small neutrino masses. It also provides a dynamical mechanism for baryon 
asymmetry generation via the lepton number violating decays of the N [4]. The induced lepton 
asymmetry is converted into baryon asymmetry via the electroweak sphalerons [5] (for reviews of 
leptogenesis see Ref. [6,7]). 

Attractive as it is, the SUSY seesaw framework is not without its problems. First, the generic 
leptogenesis mechanism is impossible to test experimentally. This is primarily because the dynam- 
ics occurs at a very high energy scale, beyond reach of foreseeable experiments. The parameters 
that are relevant for leptogenesis are not the same that appear in low energy neutrino oscillation 
experiments. (The number of low energy observables in neutrino sector is nine, while leptogeneis in 
the general setting involves a total of eighteen parameters.) Second, in supergravity models, suc- 
cessful leptogenesis is in conflict with the gravitino abundance. This is because of the lower bound 
on the lightest RHN mass M^i ^10^ GeV, (assuming hierarchical masses for N) [8] which would 
suggest rather high reheat temperature, of order 10^ GeV. This conflicts with reheat temperature 
suggested by gravitino abundance Treheat < 10^ GeV [9,10]. 

In this paper we suggest a scenario where the aforementioned problems of the SUSY seesaw 
framework are alleviated. The gravitino overproduction problem is avoided by resorting to resonant 
leptogenesis scenario [11-13] which assumes quasi-degenerate N fields. In this case the mass of the 
N fields can be as low as a TeV, consistent with successful leptogenesis, thus avoiding the gravitino 
problem. We supplement the resonant leptogenesis scenario with flavor symmetries which restrict 
the form of the neutrino Yukawa coupling matrices. Such flavor symmetries are anyway needed to 
guarantee the near degeneracy of the states. We identify three possible textures for the Dirac 
Yukawa couplings of the neutrinos that yield two quasi-degenerate A^ fields and a sum rule for the 
neutrino oscillation angle ^13. Interestingly, in all three models, there is a single phase that controls 
cosmological CP asymmetry and CP violation in neutrino oscillations. We are able to constrain 
the range of the CP violation parameter \S\ from cosmology. Somewhat similar classiflcation of 
textures has been recently pursued in Ref. [14] and earlier in Ref. [15], [16]. Our emphasis is on the 
connection between cosmological CP asymmetry and CP violation in neutrino oscillations. It turns 
out that, in our framework, there is a lower limit on the SUSY parameter tan/5 > 12. This arises 
since the mass splitting between the qusi-degenerate A^ fields is generated from renormalization 
group flow, which depends on tan (3. 
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In our analysis we use the results of a global fit to the neutrino oscilation data [17]: 
l^^atml = (2-18 - 2.64) X 10-=^eV^ (2a) , Am^^j = (7.25 - 8.11) x IQ-^eV^ (2a) , 

sin^ 023 = 0.39 - 0.63 {2a) , sin^ = 0.27 - 0.35 (2a) . (1) 

Currently we do not know the sign of Am^^^^, i.e. whether neutrinos have normal mass hierarchy or 
inverted mass hierarchy. Also, the value of the third mixing angle ^13 is unknown. Only an upper 
bound [17] 

sin^ ^13 < 0.04 (2a) (2) 

is available currently. Nothing is known about the CP violating phase S (and also about two 
'Majorana' phases) of the leptonic mixing matrix. 

We will identify explicit models wherein these unknown mixing parameters are significantly 
constrained. It will be highly desirable to relate the CP violation parameters in the leptonic mixing 
matrix with the cosmological CP asymmetry. Such a strategy was pursued successfully in Ref. [18]. 
While in Ref. [18] a close connection between cosmological CP violation and neutrino CP violation 
was realized, since the setup used hierarchical RHN masses, straightforward SUSY extension of that 
scenario would lead to gravitino overproduction. Our texture models are tailor-made for resonant 
leptogenesis, which would avoid this problem.^ 



2 Texture Zeros for Predictive Models 

Let us consider the lepton sector of MSSM augmented with two right-handed neutrinos (RHN) A^i 
and A''2. The relevant Yukawa superpotential couplings are given by 

with We^fYee'hd, W.^FY^NK-^N'^MnN, (3) 

where and hu are up and down type MSSM Higgs doublet superfields respectively. We will work 
in a basis in which the charged lepton Yukawa matrix is diagonal: 

ye = Diag(Ae, A^, A,) . (4) 

As far as the RHN mass matrix is concerned, we will assume that at high scale (identified with 
the GUT scale later on) it has the form 

^^=(1 o)^- 

This form of Mjv is crucial for our studies. It has interesting imphcations for resonant leptogenesis 
and also, as we will see below, for building predictive neutrino scenarios. Specific neutrino models 
consistent with resonant leptogenesis with a texture similar to (5) was investigated in [19]. Here we 
attempt to classify all possible scenarios with degenerate RHNs which lead to predictions consistent 

^For a concrete demonstration within predictive model see [19]. 
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with experiments. Thus, with a basis (4) and the texture (5) we can discuss possible texture zeros 
in the matrix Yi,, which is of dimension 3x2. One can easily verify that two (and more) texture 
zeros in do not lead to results compatible with the neutrino data. However, with only one texture 
zero, there are scenarios compatible with experiments and leading to interesting predictions. 

The matrix Yi, contains two columns. Since due to the form of Mn (5) there is exchange 
invariance Ni ^ N2, N2 ^ Ni, it does not matter in which column of Y^, we set one element to 
zero. We choose here the second column of Y^, having one texture zero. This leads to the three 
following possible forms for Y^: 



Texture A : = 02 62 , (6) 





Texture Bi : Y^ = \ a2 , Texture B2 : = ( ^2 ^2 ■ (7) 




A few words about the parametrization, used in (6) and (7), are in order. With the basis (4) and 
the form of M^v given in (5), the one texture zero 3x2 matrix Y^, has only one physical phase. 
Other phases can be rotated away by proper phase redefinitions of the fields. Moreover, in Y,^ there 
are five real parameters 1 01,2,3 1 and two absolute values of the 6-entries. The mass parameter M 
in (5) is in general complex, but its phase is not relevant for the physics of neutrino oscillations. 
These systems lead to predictive scenarios with texture A corresponding to normal mass hierarchy 
and textures Bi and corresponding to inverted mass hierarchy. We will study these cases in 
turn. 



2.1 Texture A: Normal Hierarchical Case 

We will discuss this case in details. With (5), (6) and using the seesaw formula for the light neutrino 
mass matrix M^, — {h^'^Y^Mj^^Yj , we arrive at 

(0 0162 01&3 

0162 20262 0263 + 0362 
0163 0263 + 0362 20363 

where w ~ 174 GeV. The matrix in (8) is rank two and leads to the one massless neutrino and two 
massive neutrinos labeled m2 and m^. This structure corresponds to the normal hierarchical case, 
i.e. 

M,d-s^Diag(0, m2, m^) , (9) 

with 777,3 ^ ?772. Prom (8) we can see that the mixings 612 and ^23 are generated. The absolute 
value of the overall factor a^h^{v sin /M determines one mass scale, say the value of 7773. Besides 
this overall factor the matrix has four parameters: one phase and three real parameters. Three of 

these parameters can be fixed from three observables 6*12, 6*23 and f "^J'"' (where Am^^j = m\ and 

^"^atm ~ ''"3 ~ ^i)- condition mi = we will still have one prediction (independent 

from the value of the phase), which determines the angle ^13. 



I (vsin, 
I M 
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One physical phase remains undetermined. Indeed this single phase will be directly related to 
the CP violation in neutrino oscillations and in leptogenesis. We will discuss this connection in 
more details in Sect. 3. 

Now, let us derive the prediction of this model. To achieve this and also get other useful relations 
we will use the equality 

M,, = PU*P'M^'''^U^P, (10) 
where U is the lepton mixing matrix, given in a standard parameterization by: 

C13C12 C13S12 

U = I -C23S12 - 5235130126"^ C23C12 " S23Sl35l2e"^ S23C13 | , (H) 
S23S12 ~ C23Si3Ci2e*^ —S23C12 — C23Si3Si2e*'^ 

with Sij — sin 9ij and = cos 9ij. The P and P' are diagonal phase matrices P — Diag (e*'^i, e''^^^ e'^^^), 
P' — Diag (1, e*''^). Phases in P can be removed by field redefinition, while P' is physical, and 
contains the two Majorana phases. The matrix equation (10) gives six relations. One of them, 
namely the relation for the (1, 1) elements of M^, and the right hand side of (10) with the form of 
U given in Eq. (11), gives 

tan ^13 ~ sin 6*124 /— ■ (12) 

V "^3 

Since this case corresponds to the normal hierarchical neutrino mass spectrum (with mi = 0), 
with tlie hell) <)f (1) we have at 2a level m2 = y^Am^ ~ (8.51 - 9.01) • 10"^ eV and 777.3 = 
•\/|AmatmP + A"^soi — (^-"^ ~ ^-2) ■ 10~^ eV. Using these values in (12), together with 2a accuracy 
value of ^12, we obtain range sin^ ^^13 ~ 0.042 — 0.062. This fits well with an upper bound, within 3a, 
given in Ref. [17], while the low limit (0.042) is pretty close to the 2a upper bound of 6*13. Future 
measurements of the ^13 will test the validity of this scenario. One more word about the neutrino 
sector: since the (1, 1) element of the light neutrino mass matrix vanishes, the neutrino-less double 
/?-decay (Oz/2/5) does not take place in this scenario. That is, 777^/3 = |f/f2^2e*^ + f^fa^^al = 0. There 
is only one Majorana phase, since mi = 0, which is p — p2 — pi- This is determined from the phase 
5 as follows 

p = 7r-25. (13) 



2.2 Textures B\ and B2 : Inverted Hierarchical Cases 

The textures Bi and B2 both lead to the inverted hierarchical neutrino mass pattern. Using these 
textures (7), the form of Mat given in (5) and the seesaw formula for the light neutrino mass 
matrices we obtain: 

2ai6i 0261 "1^3 + as^i \ / • 
For Texture Bi : = I 0261 0263 , (14) 



aih + 0361 0263 20363 



M 



2ai6i 0162 + 0261 «3&i \ c^g-j^^\5 
For Texture B2 : = I 0162 + 0361 20262 O362 \- — tj-^ ■ (15) 



0361 0362 



M 
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In order to derive predictions for both cases, we can still use the relation (10), which is general, but 
for My use the forms corresponding to the cases -81,2, and for M^'^^ an inverted hierarchical form: 

M^'^s = Diag(mi, ms, 0) . (16) 

We use the same form as before for the phase matrix P, while for the P' we use P' — Diag (e'''^ , e*^^ 1) 
For cases Bi and B2 the predictive relations emerge by equating the (2, 2) and (3, 3) elements re- 
spectively (which are zero) of the expressions at the both sides of Eq. (10). Doing so we arrive 
at: 

„ ^ 9^1 sin tan ^^23 COS 5 1 AmL , 
For texture Bi : sin^ ^12 ^ - - 1 , " . , " ^ + - . ^ f . , (17) 



2 I tan^ sin^ ^13 + e^^^ | 8 | Am^ 



atm I 



^ ^ 9^1 sin 6'i3 tan 62^ cos 5 1 Am^^, , . 

For texture B2 : sin^ ^12 ~ - + 1 tt-^ 4 ^ + - , . 9 , • (18 

2 |tan2^23 + sin2^i3e2'5| 8 lAm^^^^l ^ ' 

As we see, for both cases, the deviation of sin^ 612 from 1/2 (i.e. deviation of 6*12 from 7r/4) is due 
to the non-zero value of ^13 and it also depends on cos 5^. In fact, the product sin ^13 cos 5 should 
not be too small, otherwise the angle 9i2 will be close to 7r/4 which is excluded. Using the current 
experimental data (1) (within 2(T-deviations) we obtain the following constraints for ^13 and cos 5: 

For texture Bi : ^13 > 0.12 , cos5 > 0.573 {\5\ < 0.96) , 

For texture B2 : ^13 > 0.129 , cos5< -0.614 (|7r - 5| < 0.91) . (19) 

The last terms in Eqs. (17), (18) are practically unimportant for the neutrino sector, but as we will 
see in section 3.1 they become crucial for the leptogenesis CP violation. The leptonic asymmetry 
will be determined by a phase oc — sin 6 which would vanish in the limit j-r — |^ 0. 

By the fixed model parameters (see sect. 3.1 for relation between Yukawa couplings and the 
angles 6ij, 6) we can compute one more observable. In contrast to the normal hierarchical neutrinos 
(corresponding to the texture A), cases Bi and B2 have non-zero /5/3oi/ amplitudes, for both cases 
given by 

= |t^el"^l^*^ + ^e2"^2| , with p = ^2 " Pi ■ (20) 



For 171/313 and p for scenarios Bi and B2 respectively we derive: 
For texture Bi : m^/j = y^^Am^^^c- 



2 I 2 |tgl2 - 1 + 2tgi2tg23Si3e* 



|tgi2+tg23Sl3e*^P 

P ^ tg23(l + tgi2)gi3 sm6 

2 tgi2(l - tgiasfa) + tg23(l - tgf2)si3 COS (5 ' 

|tgl2 - 1 - 2tgi2Ctg23gl3e^'^| 
|tgl2 - Ctg23Si3e^'*|2 



For texture B2 : mpj3 = J \Ami^^\c^^ 



^Similar relation has been obtained in Ref. [19] within a specific model with 6*23 — tt/4. Here, since 623 is not 
fixed from the model, we will have somewhat wider allowed ranges for ^13 and especially for 5. Cases of texture 
zeros giving these relations have been identified recently in Ref. [14]. Correlation similar to Eqs. (17) and (18) have 
been obtained within scenarios with 'quark- lepton complementarity' [20] . 



6 



cot P- = ^ tg,3(l + tg^.K3sin^ 

2 tgi2(tg^3 - sfg) - tg23(l - tgfajsia cos 5 

where tgy = tan 9ij and ctgy = cot%. Applying allowed ranges for 6 and 6*13 given in Eq. (19) and 
the measured neutrino oscillation parameters (1) (within 2a) for m^j/s we obtain: 

For textures Bi & B2 : 0.013 eV < < 0.023 eV . (23) 

Upper bounds for mpp are obtained for \5\ — 0.96 and \'k—5\ — 0.91 for cases Bi and B2 respectively, 
while lower limits correspond to 5 — and S = dz7r. Planned experiments will certainly be able 
to test viability of these predictions. Note that the textures Bi and B2 in the neutrino sector give 
results which are practically indistinguishable (besides the allowed ranges for 6). However, as wc 
will see in the next section the scenario B2 fails to generate sufficient leptogenesis, while the texture 
Bi (and also the texture A) will work very well for this purpose. 



3 Resonant Leptogenesis 

Within the scenarios considered in the previous section, we have assumed an off-diagonal form for 
the RHN mass matrix M^. This gives the desired degeneracy between the two RHN states. The 
degeneracy will be lifted with small corrections to the (1, 1) and/or (2, 2) elements of M^. Even in 
the unbroken SUSY limit, 1-loop corrections (corresponding to the wave function renormalization) 
will split the degeneracy. The SUSY breaking effects has dramatic impact on the degeneracy of 
the scalar components of A^i,2 superfields. This is discussed separately in the Appendix. As far 
as the fermionic RHN sector is concerned, the degeneracy there holds with pretty high accuracy. 
Therefore, this is an appealing framework for resonant leptogenesis, in which enhancement of the 
CP asymmetry happens because of quasi-degenerate RHN neutrinos [11-13]. One nice property of 
the resonant leptogenesis is that, it avoids the lower bound {Mj^^ ^10^ GeV) for the lightest RHN 
mass. This bound, called as Davidson-Ibarra bound, emerges within most of the scenarios with 
hierarchical right-handed neutrinos [8]. Once this bound is avoided, the reheat temperature can be 
sufficiently low to avoid the gravitino problem, which is common for low scale SUSY models [9, 10] 
with the gravity mediated SUSY breaking. 

Since our models of neutrino masses and mixings are predictive and involve very limited number 
of parameters, we expect that we will not have much freedom in the calculation of leptogenesis. As 
we have already mentioned, an important ingredient for the resonant leptogenesis is the form of 
Mjv given in Eq. (5). Note that the mass matrix of the fermionic RHNs coincides with Mn of the 
superpotential mass term. First we will discuss radiative corrections to the superpotential mass 
matrix Mjv, which directly can be applied to the fermionic RHNs. This structure can be justified by 
some symmetry at high scale. However, at low energies, due to the radiative corrections the (1, 1) 
and (2, 2) entries in Mjv will receive non-zero corrections. These corrections are calculable thanks 
to the well defined neutrino models we have presented above. To be brief, eventually two RHNs 
are become quasi- degenerate, and the CP asymmetries ei and 62 generated by out-of-equilibrium 
decays of the fermionic components of A^i and A^2 states respectively are given by [12, 13] 



lm{Yjt)l, (M| - Mf ) Mir2 
{YjY,)n{YjK)22 {Mi - Mlf + Mfr^ 



61= V^.2-^^H,^^n^2 e2 = 6i(1^2). (24) 
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Here Mi and M2 (we will use the convention M2 > Mi) are the mass eigenvalues of the RH neutrino 
mass eigenstates. 1^ = Yi,Un is the Dirac Yukawa coupling matrix in the basis where RH neutrino 
mass matrix is diagonal and real: UJ^MnUn — Diag (Mi, M2). is the tree level decay width of 
Ni (mass eigenstates of RHN) and is given by Fj = ^{YjY^)ii. From (24) we see that in order to 
have non-zero CP asymmetry two conditions need to be satisfied. First, the RHN masses should be 
split, and secondly the element {YJYj,)i2 must be complex. To realize both of these conditions, we 
need to include radiative corrections into our study. As we will see shortly, the desired result can be 
obtained only at two-loop level. In our treatment we assume that the textures we have considered 
are realized at the GUT scale Mq ^ 2 • 10^^ GeV. At low scales, due to the renormalization group 
effects the zero entries in the flavor matrices will receive some corrections. To compute these 
corrections we set up the RG equation for the matrix Mjv (only its renormalization has relevance 
for us), which at two-loop order is given by [21]: 

IGtt^^Mjv = 2MnY^Y, + 2YJy:Mn 
\Mn {YjY.Y^Y, + Y^YX + YX{^>^t + tr{YjY^)) 



87r2 



^ (If y^y/y; + yjy;yjy; + yJy;(3A? + tr{YX)) + ^^m^, {gi + bgi) , (25) 



g^2 x v e e v ' v v v v ' v v \ t ' \ v >^ ) ) "^jv i 2077^ ' 

where i = In The first line in (25) corresponds to the 1-loop correction and will be responsible 
for the mass splitting between RHNs. However, the two-loop correction, presented in a second hne 
of Eq. (25), will be crucial for the CP phase of {Y},Yi^\2- Since we intend to have Mi^2 ~ 10^ GeV, 
in order to get reasonable scale for the light neutrino mass, the matrix elements of Yy should be 
much less than unity. Thus, we can solve the RG equation analytically to a good approximation. 
One-loop correction to the Mat can be found from (25) to be 

5Ml,-'°- ^ {M^yX + 1- ^ • (26) 

From this we see that at scale /i = M the form of Mjv will become 

Interestingly enough, this structure, of correlated phases of (1, 1) and (2,2) entries of M^v, persists 
also at two-loop order. What is more important, one can sec that at the one-loop level the phase 
of is determined by the phase of {YJY^) 12 and therefore (Y^Y^) 12 will be real at this level. This 
property can be easily seen also from different angle. Regardless of the form of Y^, (including all 
possible radiative corrections to it), it can be written in the form 



Yy^U \ a2 P , (28) 




with U some unitary matrix, 01,2,^3 all real parameters and P — Diag(l,e*^). Using now the 
form (28) in the first line of (25), one can show that U drops out and we remain with the non 
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physical phase ^ which can be absorbed in With this, any complexity in 5M^~'°°^ and (YjYi,)i2 
disappears and we have no CP violation at the one-loop level. That is why it is important to include 
two-loop radiative corrections for the renormahzation of M^. Indeed, the term MnYJY^Y^Yi, in 
the second line of Eq. (25) is important. The appearance of the combination Y^Y^ plays an 
important role. With the basis (28) wc sec that the matrix U does not disappear, and thus we 
expect to have CP violation (induced through two-loop correction). From (25), this correction can 
be approximated as follows: 

SMl'''"'' ^ {M^YXY^ + YJy:yJy:M^) ^^^^ i?. In ^ . (29) 

where we have suitably absorbed CP conserving and flavor universal corrections (coming from the 
entries Tr(yjy^), gf, \^ etc.) in the overall scale M. The RG factor {ii = {e,fj,,T)) is for the 
running of Yg Yukawa couplings and can strongly deviate from one. It is defined as: 

^'"^"A|,,(M.)ln^- ^'"^ 

In the approximation (29), the fourth powers of Y^, have been neglected. Actually, for calculating 
the mass sphtting M| — Mf - the combination appearing in (24) - it is enough to keep only correction 
5M^~'°°'' of (26). However, to deal with the CP violating effect, we need to include also two-loop 
effects. Thus, at the scale /i — M ior we use 

Mr,{M) = Mjv(Mg) + 5M^-'°°P + 5M^-'°°p , (31) 

with M]\i{Mq), 5Mj1^~'°°'' and 5Af^~'°°^ given by Eqs. (5), (26) and (29) respectively. This completes 
the calculation of supersymmetric part, which will be useful for calculation of leptogenesis via 
fermionic RHN decays. However, inclusion of soft SUSY breaking terms, in general, may affect the 
leptogenesis induced through the right-handed sneutrino decays. In an Appendix we studied this 
case separately and shown that under plausible assumptions the right-handed sneutrino decays 
practically do not contribute to the net baryon asymmetry. Thus, we should relay on the fermionic 
RHN decays which, as we show below, generate sufficient baryon asymmetry. 

3.1 Asymmetry Via Fermionic RHN Decays 

Leptogenesis for Normal Hierarchical Case 

For this case we will take the form of Y^, given by Eq. (6). For leptogenesis study, it is convenient 
to parameterize this Yukawa matrix as follows: 



Texture A : Y^ = xa2 h , (32) 




where the couplings q;i,2,&, /? and x are real parameters. Only single phase (p appears. This has 
been achieved by suitable redefinition of phases of Zi,2,3 and Ni^2 superfields. First we will relate the 
parameters appearing in Y^ to some observables. The relation (10) enables us to express q;i,q;2,6 
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and j3 in terms of x, neutrino mass, the scale \M\ and lepton mixing matrix elements. Also (j) can 
be determined by the phase 5 and the leptonic mixing angles. Doing so, we arrive to the following 
relations 






s 



Figure 1: Correlation between and 5. Left side: normal hierarchical case (texture A). Right 
side: inverted hierarchical case (texture Bi). The vertical lines, for right panel, correspond to the 
maximal allowed value of \5\ = 0.96. 



U^2 ^e2 Ue2 'J e2 

These will be useful upon studying the leptogenesis. As we have already mentioned, remarkable 
thing is the fact that there is a single CP violating phase (f) which is related to the phase 5 controlling 
the CP violation in the neutrino oscillations. The same phase will appear in the CP asymmetry of 
the resonant leptogenesis. In Fig. 1 we show correlation between and 5. 

Furthermore, applying expressions (26), (29), for the splitting parameter 5n of (27) we obtain 



\ ba2 + e^^ -^eAf,\n^^ , (36) 



where we have ignored the couplings Ae and because the main effect is obtained by the tau 
Yukawa coupling. In Eq. (36), the coupling A^ is defined at Mz scale, and therefore the quantity 
Rt- accounts for the renormalization effects mostly due to A,- running, and is given in Eq. (30). Now 

-Di, 
N 



we can give the unitary matrix Un diagonalizing Mjv (by the transformation UJ^MnUn = M^^^* 



^iv ^ ^ ) ' (37) 
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X 

Figure 2: Baryon asymmetry for normal hierarchical case (texture A), for different values of tan/3 
and M = 10^ GeV, 5 = 1.3. 



where 



rf = ATg(ba, + e'^-^e''f'] . (3^ 



At the same time we have 

{YJY^)2i = P^x \ba2 + e^'^l e*"' , with r]' = Arg {ba2 + e'^) . (39) 

Therefore, the complex phase appearing in iYjY^)2i will be proportional to the mismatch t] — r]', 
which using (38) and (39) takes the form 

, XlRr ba2 . , 

Note once again that the phase t] — rj', determining the lepton asymmetry, is proportional to sin 0, 
which itself is related to the phase 6 of the lepton mixing matrix. The model gives the relation 
between them by Eq. (34). Also, it is rather impressive that other parameters, b and a2, appearing 
in (40) can be calculated by the lepton mixing matrix elements through the relations (33), (35). 
The masses of two right handed neutrinos are 

•A , B"^ Mr 

Mi = \M\{1-k), M2 = \M\{1 + k), with n:^\x{ba2 + e"^)\j-^\n-^. (41) 

Here, the unknown parameter x appears which is free and can be varied upon numerical calcula- 
tions. Finally, we give expressions build from the elements of the matrix (YJYu) appearing in the 
expressions of the CP asymmetries of Eq. (24). These are: 

iYj%)u ^ Y (x'(l + a? + al) + 1 + + 2x\ba2 + e^l) , 
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(yjy,)22 ^ y + a? + a^) + 1 + 6' - 2x|&a;2 + e^"^]) , 

Im(F;n)^, c (^'(1 + + «D - 1 - ^^^f^ ■ (42) 

In order to compute generated baryon asymmetry of the Universe, recall that the Icpton asymmetry 
is converted to the baryon asymmetry via sphaleron processes [5] and is given by ^ ^ —1.48 x 
10~^(fi;/^^''ei + fi;/*^^^e2), where the efficiency factors Kf^^'"^^ are given by the extrapolating formula [6]: 

(1,2) ^ ( 3.3 X 10-=^eV / mi, 2 ' ''' ' 

'^^ y mi,2 VO-55 X 10-3eV 

with™. = (li^(ftyj„, ft, . . (43) 

Collecting all this, we can now calculate — . One can try the different values of M in a mass 
range which would not cause the gravitino problem. We can also try different values of the phase S, 
relevant also for the CP violation into the neutrino oscillations. As we have already mentioned, there 
is one more free parameter x, which we will vary. It is quite interesting that this system, by requiring 
to have baryon asymmetry in the range of the observed amount — (8.75 ± 0.23) • 10~^^, 

dictates the preferred range for the MSSM parameter tan/5. The reason for this is simple. The 
strength of the Yukawa coupling A,-, determining the amount of the CP violation [see Eq. (40)], 
depends on the value of tan/?: A,- = — a/1 + tan^ p. By simple but quite complete numerical 
simulation we obtain, in this model, the low bound on the tan/3. Upon calculations we take into 
account the renormalization effects. Namely, the running of A,-. Obtained low bound for tan/3 is: 
tan/3 > 12 (corresponds to \6\ ~ 1.3 and M = 10^ GeV, Rr = 0.617). Smaller values of tan/3 
do not give sufficient baryon asymmetry. This also indicates that the non SUSY version (i.e. SM 
augmented by two RHNs) of this scenario will fail to generate baryon asymmetry through the 
leptogenesis. The presented scenario also allows to derive the low bound for the absolute value of 
the phase 5. This comes out from the maximal allowed value of tan /3 < 58 (from the requirement 
that Xb^r ^ 1 all the way up to the GUT scale). With tan/3 = 58, M = 10^ GeV {Rr = 2.17) 
in order to have needed baryon asymmetry we should have \d\ > 0.012. It is interesting to note 
that for tan/3 ^ 35, for generating the baryon asymmetry we need \6\ ^ 0.1. This limit for the CP 
violating phase is within the reach of future experiments. In Figs. 2 and 3 we plot ^ for different 
choices of the model parameters. 

Leptogenesis for Inverted Hierarchical Case 

Now we study the leptogenesis for the inverted hierarchical case. We note right away that the 
scenario with texture B2 of (7) does not work for the leptogenesis. The reason is following. Due to 
the zero in (3, 2) entry of this texture, it is easy to see from Eq. (29) that the A,- couphng do not 
contribute to the CP asymmetry induced at 2-loop level. The couplings Ag and A^^ do contribute, 
however they are small and can not induce needed asymmetry. 

Thus, we focus here only on case with texture Bi. For this case, it is convenient to write 
with the parameterization 

Texture Bi : ^ \ xa2 \ -(3 , (44) 
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where, as in case of texture A, by suitable phase redefinition of /i,2,3, -^1,2 superfields we left with 
only single phase 0. Remaining parameters are real. First we will express the model parameters 
<^i,2,b,(3 in terms of matrix elements of U, neutrino mass m2, the M, and x. By solving the 
equations derived from the relation (10) we obtain 



For Texture Bi 



B1B3 



a2 



Bo 



Bi 



B. 



B2 



m2 



BiM 



X 



Arg 



BIB4 

bM 



where 
Bi = 



iP - 



'^rl rj2 



B2 



Uf^2Ur2~Url 



B. 



Ue2Uii2~Uel 



'-'fj.2 



f/, 



B, = Ui,-Ut, 



2 ^ 



1/2 



(45) 



(46) 



As we see, also in this case the phase (f) is related to the 5-phase of the leptonic mixing matrix U 
(see Eq. (11)). In particular using the relation (17) in (45) for cf) and performing simple algebra we 
derive 



Arg 



lAmftrr,! 4si3 tan 6*' 



- 1 



723 



cot 



23 



Amft^l 4sin0 



sin 5 
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(47) 



Since the phase will appear in the leptonic CP asymmetry, with relation (47) we will be able to 
make calculations in terms of measured neutrino oscillation parameters and the CP phase 5. In 
Fig. 1 the correlation between and 5 is shown. 

Now we are ready to investigate the leptogenesis for the inverted hierarchical scenario (-Bi). The 
way of calculation is same as was presented in the previous subsection, so we will keep discussion 
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short and give only several expressions and final results. Using the expressions of Eqs. (26)-(31) 
and the form of the texture Bi in (44), for the phase mismatch (in analogy of Eq. (40) we obtain 

~ sine/. f48) 

where here and below we will use superscript '(-Bi)' in order to distinguish expressions corresponding 
to the scenario Bi from those of the texture A. Moreover, for the splitting parameter (in analog to 
Eq. (41)) we have 

iQ2 T\/r 

K^^^) c^\x{hai + e'^)\^hi-^ . (49) 

47^^ M 

We will also give the expression for \miY^Yy)\^ which will help to understand some physics. We 
have 

(lm(y<K)L)""' ^ (-'(1 + + a|) - 1 - f^) 1^;^^ ■ (50) 

Note, that according to (47) the phase is close to vr and one may suspect that also final result 
for the CP violation should be suppressed by the factor~ AmgQi/(|Am^^jjj|4sin^i3) f» 1/20. How- 
ever, such suppression do not takes place because the combination \hai + e*"^!, appearing in the 
denominator of the last multiplier of (50), is suppressed by precisely same factor! Indeed, using the 
relations of Eqs. (17), (45)-(47) we derive 

|M,+e-|=.^i2^. (51) 

With these for the combination appearing in (50) we get 

• smd, (52) 



\bai + e'^\ 
14 




showing that suppression factors mentioned above drop out and it is maximized with \6\ 1.115 
(maximal allowed value which is acceptable for viable neutrino sector). Moreover, because of 
the suppression of the combination |6a;i + 6**^1, also the RHN mass splitting parameter in (49) gets 
additional suppression, which makes two RHNs more degenerate. This also gives some enhancement 
of the (resonant) CP asymmetry factors ei,2- 

Without bothering to give other expressions, we will move to the presentation of the main 
results. In this scenario, from the requirement of needed baryon asymmetry, the tan f3 is bounded 
from below. Interesting thing is that the leptogenesis dictates tan/3 > 21 (lower values do not 
give sufficient baryon asymmetry) For obtaining this low bound we have taken M = 10^ GeV 
[Rj. ~ 0.71) and maximal allowed value for the 6 ~ 0.96. Note that within this scenario low values 
of M give larger lepton asymmetries. It is also possible to derive low bound for \6\. This is obtained 
by largest allowed (from requirement A;,,- < 1 up to the GUT scale) value of tan/5. Namely, with 
tan/3 ~ 58, M = 10^ GeV (we have for these choices Rr — 1.95), needed baryon asymmetry can be 
generated with \6\ > 0.021 (note that for \6\ = 0.021, for acceptable solar mixing angle we should 
choose sin^^23 — 0.6 and sin^^is ~ 0.04). Worthwhile for noting that for tan/? < 45 for generating 
sufficient baryon asymmetry we need \6\ ^ 0.1. The latter value is within the reach of planned 
experiments. We have performed numerical calculations without approximations and made sure 
that our analytical expressions, presented above, are good approximations. In Figs. 4 and 5 we 
show baryon asymmetries for several different choices of the model parameters. 



4 Summary 

In this paper we have considered an extension of MSSM with two quasi-degenerate right-handed 
neutrinos. Our motivation was to realize resonant leptogenesis which avoids the gravitino problem 
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generic for low scale SUSY scenarios. With this setup we have classified all viable texture zeros of 
the neutrino Dirac Yukawa matrices which lead to consitent predictions. We find three predictive 
scenarios, each with one texture zero. One model has normal hierarchical neutrino mass spectrum, 
while the remaining two have inverted hierarchical mass pattern. The predictive power of these 
models show up also in the resonant leptogenesis. The model with the normal mass hierarchy 
(texture A) and one of the inverted hierarchical scenario (with texture Bi) lead to the successful 
leptogenesis. In Appendix we have discussed the impact of the soft SUSY breaking terms on the CP 
asymmetry generated by RH sneutrino decays and concluded that with natural choice of the soft 
SUSY breaking terms, scalar RH neutrinos do not contribute sizably to the total baryon asymmetry. 
Thus, the baryon asymmetry is due to fermionic RHN decays and the leptonic CP phase is directly 
related to to CP violation in neutrino oscillation. Putting together the predictions from the neutrino 
sector and the results from leptogenesis calculations, we have obtained the following predictions: 

For normal hierarchical case (texture A) 

sin^ ^13 > 0.05 , \S\ > 0.012 , m^j/? = , tan/3 > 12 ; 

with tan/5<35, \6\ > OA . 

For the inverted hierarchical case corresponding to texture Bi: 

^13 > 0.12, 0.021 < < 0.96 , 0.013 eV < m^j^? < 0.023 eV , tan/3 > 21; 

with tan/3<45, \S\ > OA . 

The texture B2 do not generate the baryon asymmetry within this scenario and other mechanism 
need to be invoked [23] . However, from the viewpoint of the neutrino sector the texture B2 is viable 
and gives: 

^13 > 0.129 , \n-S\< 0.91 , 0.013 eV < m^^ < 0.023 eV . 
Future experiments will examine the viability of these scenarios. 
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Appendix: Asymmetry Via N Decays 

In this appendix we will discuss the contribution to the net baryon asymmetry from the out of 
equilibrium resonant decays of the right handed sneutrinos (RHS). With inclusion of the soft 
SUSY breaking terms, the RHS mass spectrum and couplings will be altered and one should 
expect result different from that corresponding to the fermionic RHN decays. Besides soft SUSY 
breaking couplings, there are other particularities, highlighted below, which distinguish cases of 
RHN and RHS decays. We are considering the system with two RHN superfields Ni^2 which have 
two complex scalar components iVi,2- With SUSY breaking term, the masses of RHS's will differ 
from their fermionic partners' masses. Thus we will have four real mass-eigenstate RHS's ^1=1,2,3,4 
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with masses Mj=i^2,3,4 respectively. Assuming that the SUSY scale is smaller (at least by factor 
of 10) than the scale M (the overall tree level mass for the RHN superfields) we expect that 
the states fii remain quasi-degenerate. To study the resonant n-decays we will apply ressumed 
effective amplitude technic [12]. An effective amplitudes for the real fii decay, say into the lepton 
la {a. = 1, 2, 3 is a generation index) and antilepton respectively are given by [12] 

A _ ^ _ V- %(M,)(1-^,,-) ^ _ o* _ c* n,,(M,)(l-^,,) 

where Sai is a tree level amplitude and Hj^ is a two point Green function's (polarization operator 
of fii — fij) absorptive part. The CP asymmetry is then given by 

= 7-. ^4 • (54) 



We will apply (53) and (54) for our scenario, however, also derive general expressions applicable 
for different models. 

Toegether with superpotential couplings (3) we include the following soft SUSY breaking terms 

V^^ = lA^Nhu - -N^BnN + h.c. (55) 

We do not display here soft mass^ terms, such as ?TJi,2l^i,2p, because plays much more signif- 
icant role in the sphtting of RHS masses. For simplicity we will assume at GUT scale {Ma) the 
'proportionality' A^, oc Y^, and degeneracy in B^ oc Mjv. Thus, 

aX 11 = Mq ■ Ay = niAYy , B^ = ttibMn^ . (56) 

Similarly, for the charged lepton sector we can assume Ag = m^Ye. Performing RG studies, similar 
way as we have done in section 3, we will have 

at /. = M: Bnc^ttibMI^ ^ l"'' _(1 + 2^)5^ j " ^^^^ 

Note that with A^, = rriAYi, and Ag = msYe at high scale, the Ai, will remain well aligned with Yj^ 
also at low scales. With diagonalization of total mass matrix of the RHS's, for mass-eigenstate (hi) 
masses we get 

= |M|2(1 - \dN\f - \M\\mB - {niB + 2mA)\SN\ \ , 
Ml = |M|^(1 - \5N\f + \M\\mB - (mB + 2mA)|(^iv|| , 
M| = |M|2(1 + I^Tvl)^ - \M\\mB + {ruB + 2mA)|(5jv|| , 

Ml = |M|2(1 + \5^\f + \M\\mB + {niB + 2mA)\5N\\ ■ (58) 
Interaction of h states with leptons and sleptons has the form 

hulYpfi + hJYBn + h.c. (59) 
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where 

Yf = Y,V , Yb = Y,M*^V* + A,V , where V = Ur, (p„e'^ Pd 

With these we can calculate the absorptive part of the polarization diagram with external legs fii 
and Uj. At 1-loop level it is given by 

n.,(p) = ^ {p'yIyf+p^y^y* + Y^Yb + .. , (61) 

where p denotes external momentum in the diagram. 

Now we are ready to calculate the lepton asymmetry. Note that in unbroken SUSY limit, 
neglecting finite temperature effects (T 0), the decay does not produce lepton asymmetry. 
The reason for this is following. The decay of in two fermion is A^ — > Ihu, while in two scalars 
is A" — > r/i*. Since the rates of these processes are same due to SUSY (at T — 0), the lepton 
asymmetries created from these decays cancel each other. However, with T 7^ the cancelation is 
partial and one has 

g. = eiihi lhu)^BF , (62) 

with temperature dependent factor A^^ given in [22]. We note that Eq. (62) is valid when we 
have the alignment A^, = rriAYu. Without this alignment other terms in r.h.s of (62) proportional 
to uia/M will appear. Since we are assuming the alignment and m^/M < 0.1, the SUSY breaking 
effects would not affect decay amplitudes significantly and we can apply (62) for our study. Thus, 
we just need to compute ei{hi Ihu) - the asymmetry created by fii decays in two fermions. Using 
in (53) Sai — {YF)ai, with (54) after straightforward calculation we obtain 

r ll^ 1 (Mf - M|)Im(n,,) ^ 

eAui Ihu) — — J S 2 > — lm{YJrYF)ii + 

{yIyf)u 1 Y (M^ - M/)2 + |n,,f ^ ^ 

Y-T (V^V^ (MJ - )Im(n,,) - (M| - )Im(n,,; ^ . 
} lm{Y^YF)kj . ^ — : z ^UjiUki } . (63) 

tt {{Mi - Mf)^ + |n,,f) [{Mi - Ml)- + in,. f 



^kk 



In (63) for the absorptive part IT we should use (61) with p — M^. Now, the baryon asymmetry 
created from the lepton asymmetry due to n decays is: 

4 _ 4 
^ ~ -8.46 • 10-^ V -^Vi = -8.46 • 10"^ V e,(n, ^ lhu)Vi , (64) 

where we have taken into account that an effective number of degrees of freedom, including two RHN 
superfields, is (?* = 228.75. rji are an efficiency factors which depend on fhi ~ ("'^™/^) 2(YpYf)u, 
and take into account temperature effects by integrating the Boltzmann equations [22]. Before 
discussing this in more details, it is more instructive to see what are the effects of the soft SUSY 
breaking terms in the CP asymmetry given by Eq. (63). The parameter is controlled by the 
imaginary parts of the elements of the matrix Y^Yf- First note that the phase 9 appearing in this 
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matrix (sec Eq. (60)) for M < 10'' GeV is ^ ^ 10 and can be safely ignored. With this, the 
matrix YpYp has the form 

where Y^ = Y^Un is the same matrix appearing in the CP asymmetries (24) induced by fermionic 
RHN decays. Note that the matrix a has purely imaginary entries and they can be new sources 
for the CP violation. For instance, the element {YpYp)i2 has the large phase. This means that 
there happens the 'conversion' between ni and n2 states. On the other hand, from Eq. (58) one 
can see that the degeneracy of Mf and M| is split by the i?-tcrm and unless ~ 10 MeV the 
resonant enhancement does not happen (similar to the case of soft leptogenesis [22]). Since the 
natural value of is from few x 100 GeV to few TeV, we conclude that this channel does not give 
important contribution to the CP asymmetry. For those states amongst which degeneracy is not 
ruined (the 'pairs' fii — and n2 — n^j by the B-terms, the CP asymmetry is controlled not by 
imaginary components of a but by lm{YjY^)i2 (like to those corresponding to the fermionic RHN 
decays, Eq. (24)). Thus, the CP asymmetry via decays would not be larger than asymmetry 
generated due to their fermionic partners. Moreover, due to the efficiency factors r)i, the ub/s 
turns out to get additional suppression in comparison to the ub/s (the total baryon asymmetry 
due fermionic RHNs). We have checked this on two examples corresponding to the textures of A 
and Bi. Namely, we have performed calculations for (m^, itlb) = {lO^i, 10^) GeV and for several 
choice of model parameters (tan/5,M, 5). For a given set of these parameters, for fixed x we can 
calculate the values of the masses rhi — 2{Y^Yp)ii. With given values of Wj, according to 

Ref. [22] we picked up the corresponding values of rji and with help of Eqs. (63), (64) calculated 
ub/s. For the texture A we obtained as- < 4 ■ 10^''', while for the texture Bi. ^ < 10~^. These 

' ' ^ riB 

confirm that the baryon asymmetry via n decays is a negligible effect. For completeness we also 
examined the case corresponding to texture B2. The latter does not give relevant asymmetry also 
through n decays. 



, with a — - 
2 



1 i 
-i 1 
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